Effect of electron correlations on (001) Fe/MgO interfaces. 
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We developed a parametrization of transmission probability that reliably captures essential ele- 
ments of the tunneling process in magnetic tunnel junctions. The electronic structure of Fe/MgO 
system is calculated within the quasiparticle self-consistent GW approximation and used to evaluate 
transmission probability across (001) Fe/MgO interface. The transmission has a peak at +0.12 V, 
in excellent agreement with recent differential conductance measurements for electrodes with an- 
tiparallel spin. These findings confirm that the observed current-voltage characteristics are intrinsic 
to well defined (001) Fe/MgO interfaces, in contrast to previous predictions based on the local 
spin-density approximation, and also that many-body effects are important to realistically describe 
electron transport across well defined metal-insulator interfaces. 
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Energy dissipation in switching a logic unit (transistor) 
is perhaps the most important bottleneck to continued 
realization of Moore's law scaling in the density of inte- 
grated circuits P, 0. New materials with special, well 
denned interfaces J3| and functional properties Q offer 
promising routes to circumvent losses. Most new schemes 
exploit tunneling phenomena across metal-insulators or 
semiconductor interfaces, so they will play an increas- 
ingly important role as we explore fundamental limits of 
density (miniaturization) in integrated circuits [2j. 

In thin FM/insulator/FM heterostructures, where an 
FM is a ferromagnet, spin polarized electron tunneling is 
observed. In such "magnetic tunnel junctions" (MTJs) 
the tunneling resistance changes when alignment of the 
two FM electrodes are switched from an anti-parallel con- 
figuration (APC) to a parallel configuration (PC). This 
property is encapsulated in the tunneling magnetoresis- 
tance (TMR), T MR = (G P - G a )/Ga. G p and G A are 
conductivities in the PC and APC 0, H- These het- 
erostructures are of particular technological interest be- 
cause of the recent discovery of very lar ge T MR in highly 
crystalline (001) Fe/MgO/Fe MTJs 0,0. Equally im- 
portant to the large TMR is the record low critical cur- 
rent needed to manipulate (switch) magnetization with- 
out external magnetic field, relative ease of fabrication 
and reproducibility [llj]. At a fundamental level, the 
(001) Fe/MgO/Fe system presents an excellent opportu- 
nity to investigate the role of the interface and test our 
ability to predict critical properties controlling transport 
in a heterostructure with well denned interfaces [§] . 

The local spin density approximation (LSDA) has been 
used to predict TMR in Fe/MgO/Fe system p~3l4T7j . Re- 
ports are largely in agreement with each other but at 
variance with differential conductance measurements for 
thin MTJs [H, El- The LSDA predicts a narrow band 
of interface resonance states (IRS) of minority electrons 
(Fig. 1) which overlaps the Fermi level Ep. This pre- 
diction has two consequences: first, a sharp reduction 



in TMR at voltages on the order of ^0.02 V, due to 
sharp (resonant) reduction in minority spin contribution 
to differential conductance in PC [TBI, [l8j]. It is called the 
'zero-bias anomaly'. Secondly, it predicts a rise in TMR 
at larger voltage, ~0.1 V, due to reduction in differen- 
tial conductance in APC [Hjj. Neither effect has been 
observed experimentally; instead differential conductance 
in APC increases monotonically, with d 2 I/dV 2 reaching 
a peak near 0.12 V |12] or 0.15 V Zermatten et al. 
attributed this effect to interface states, observed at the 
(001) Fe free surface by STM measurements. 

The absence of 'zero-bias anomaly' in PC is usually ex- 
plained as a slight asymmetry of the electrodes from e.g., 
disorder, which breaks matching of the resonant states. 
Thus for any applied bias the minority-spin contribution 
to current in PC is much smaller than the majority-spin 



contribution for MgO thicknesses larger than lnm [14 1. 
As regards the APC, whether inconsistencies between the 
LSDA and experiment are due to some extrinsic phe- 
nomenon (e.g. disorder), or a failure in the LSDA, has 
remained an open question. The LSDA underestimates 
bandgaps, which strongly affects the band structure at 
imaginary k governing evanescent decay. It also poorly 
describes the Schottky barrier height [23(. Thus in the 
LSDA, most of the key parameters responsible for TMR 
are somewhat suspect. 

Here we investigate electronic structure of (001) 
Fe/MgO/Fe using the Quasiparticle Self-consistent GW 
(QS GW) approximation for the electronic structure, 
which does not rely on the LSDA. We find significant 
corrections to energy and fc-space character of the mi- 
nority spin channel interface states. We investigate how 
these electronic structure corrections alter I-V character- 
istics in the APC, and account for available experimental 
results 
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QSGW uses self-consistency to minimize 
the many-body part of the hamiltonian, thus allowing 
accurate determination of quasiparticle (QP) levels with 
low-order diagrams. It has been tested for a wide variety 
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of bulk material systems and has been shown to be a good 
predictor of materials properties for many classes of com- 
pounds composed of elements throughout the periodic ta- 
ble 13" 21 1 . It vastly improves on the accuracy of any ex- 
isting DFT method, as well as standard implementations 
of GW, namely (1-shot) perturbations around the LSDA. 
Also it surmounts problems of self-consistency inherent 
in a true self-consistent GW scheme. Self-consistency is 
particularly important in ionic compounds such as MgO 
2o| . Moreover, it may be necessary for reliable descrip- 
tion of the metal/insulator interface since screening at 
the interface can be particularly important. For example, 
screening has been shown to strongly affect the molecular 
levels of benzene near a graphite interface [22[ , a correla- 
tion effect not captured by a Kohn-Sham theory. As we 
show here, QSGW applied to this highly heterogeneous 
system appears to have the same uniform accuracy found 
in homogeneous materials systems. 

We develop a formula to obtain the transmission prob- 
ability T that avoids direct calculation of the tunneling 
via Landauer-Buttiker theory. T is parameterized by the 
local density of states (DOS) inside the tunneling layer. 
As we will show, it does an excellent job at reproduc- 
ing the full Landauer-Buttiker transmission for a given 
one-body hamiltonian; it thus makes possible predictions 
of transport within the QS GW approximation. Used in 
conjunction with QS GW, T calculated in APC is in ex- 
cellent agreement with observed I-V characteristics. This 
demonstrates that the observed TMR is not an artifact of 
imperfections at the interface, but an intrinsic property of 
it. An important corollary is that many-body corrections 
to the LSDA for electronic structure can have profound 
effect on transport properties. In this particular case the 
dominant correction to the LSDA is a shift in the IRS, of 
the same magnitude as a typical bias voltage. But gener- 
ally speaking we can expect QS GW to describe electronic 
structure in inhomogeneous systems with vastly better 
accuracy than commonly adopted approaches and thus 
investigate implication of these corrections for transport. 

Owing to heavy computational costs, the Fe/MgO in- 
terface is modeled in QSGW with a periodic slab of 5 
Fe and 5 MgO layers ordered on the (001) plane. All re- 
sults reported here adopt a gene ralized linear muffin-tin 
orbitals (LMTO) method [2JJ,[2l|, and use the relaxed nu- 
clear coordinates of the Fe/MgO interface from Ref. [24j ]. 

Minority DOS of the Fe/MgO superlattice, projected 
on the Fe interfacial atom, and onto the oxygen atom 
in the second layer from the interface, are depicted in 
Fig. Q] For both LSDA and QS GW two narrow peaks 
are seen, separated by ~0.15 eV. The iron and oxygen 
projected DOS have similar shapes, indicating that these 
peaks originate from the same interface resonance states. 
These states fall near midgap in MgO and decay expo- 
nentially, the decay being more pronounced in QSGW. 

The LSDA puts one peak just at Ep, in agreement 
with several prior LSDA calculations |13l4l7j ). There is 




FIG. 1: (color online). Minority electron DOS projected to 
the surface Fe layer (left panel) and the oxygen atom in the 
second layer from the Fe/MgO interface (right panel). QS GW 
results are depicted by a (red) solid line, LSDA by a (blue) 
dashed line. Note that the panels have different scales. Black 
dotted line dep icts \d 2 I/dV 2 \ measured in the antiparallel con- 
figuration [13 ] (arbitrary units). 



a corresponding IRS resonance in QSGW, but it falls at 
-Ef+0.12 eV. This difference is important: it is compara- 
ble to typical bias voltages. 

Through calculation of T we can establish a connec- 
tion with I-V measurements (black dotted line in Fig. [I). 
The observed \d 2 I/dV 2 \ has peaks at both +0.12 V and 
—0.12 V, with the former being more pronounced. The 
difference can be explained by the structural asymmetry 
[l2l ]: the top interface is grown last and is rougher than 
the bottom one. In the following we show quantitatively 
that a combination of featureless DOS in the majority 
channel (not shown) and sharply peaked IRS states in 
the minority channel explains the observed peak in the 
differential conductance \d 2 I/dV 2 \ at +0.12 V. (Voltage 
is defined so that forward bias samples minority states of 
the bottom electrode with E>Ep). We use the following 
expression for T at zero applied bias: 

T°,{E) = A^DOS£ B (i) e - 2 ^ d -DOS^(i?) (1) 

N\ 1 is the number of k-points in the 2D Brillouin zone nor- 
mal to the interface, a and a' denote spin polarizations 
of the left and right electrodes, DOS£ s (L) is the DOS of 
electrons with spin a projected onto nucleus L, chosen at 
will somewhere in the MgO close to the left Fe/MgO in- 
terface, and DOS^^i?) is the corresponding DOS of spin 
a', projected onto a nucleus R in MgO close to the right 
Fe/MgO interface. 7k£ is the smallest (spin-independent) 
imaginary wave number for evanescent states inside the 
MgO barrier with given k and energy E—& property of the 
complex band structure of bulk MgO. cIlr is the spacing 
between planes containing atoms L and R. Eqn. ([IJ ne- 
glects parallel channels with larger imaginary wave num- 
ber, which is always satisfied if the barrier is thick enough. 
It becomes exact for one-dimensional case if the projected 
DOS is replaced by the square of the wave function tjj 
propagating in correspon ding electrode that is normal- 



ized to carry unit flux 28, 29]. Since we replace the flux- 
normalized |V>| 2 with a local DOS, we include factor A to 
correct for the (effectively unnormalized) ip. Once local 
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DOS and 7ks are given, T can be calculated for arbitrary 
MgO thickness. Though Eq. ([1]) bears a superficial resem- 
blance to Julliere's formula, the latter takes into account 
only the spin polarization of electrodes, while Eq. (jTJ) ac- 
counts for barrier-electrode coupling and evanescent de- 
cay as well. These contributions are the essential ones in 
the Fe/MgO system. 

To evaluate the trustworthiness of Eq. ([TJ, we compare 
it to a complete calculation of transmission within the 
Landauer-Buttiker formalism for a system with 4 MgO 
layers and semi-infinite Fe electrodes. For this purpose we 
use an implementation within the tight-binding LMTO 
method and Atomic Spheres Approximation (TB-LMTO- 
ASA) [25LI26I. [As we show below the ASA is reasonably 
close to, but not identical with the full potential (FP) 
LSD A result; it matters little here since our purpose is to 
evaluate the reliability of Eq. (jTJ)] _ 

We first calculate transmission in APC, Tap(E), 
within the Landauer-Buttiker formalism, and compare to 
T^ P (E) from Eq. (TTJ, staying within the TB-LMTO-ASA 
method. For T® P (E) we employ DOS for a barrier con- 
taining 12 MgO layers (that ensures that local DOS is 
converged with thickness), and use for L and R respec- 
tively the O atom in the third layer from the left inter- 
face (£=03) and the second layer from the right inter- 
face (R=Of). Thus the two DOS entering into Eq. (TTJ 
are DOS^(O^) and DOS^ n (Of). L=Oj and #=Of 
because oxygen atoms have more valence electrons and 
are located closer to the interfacial Fe atoms; they thus 
better represent electrode-barrier coupling than do the 
Mg atoms. Also, for the 4-layer MgO barrier L=0^ and 
R=C>2 is the same atom. Then dpp=0 and factor j^p 
in Eq. (TTJ is not needed. As Fig. Ufa) shows, T% P (E) is 
nearly identical with the Landauer-Buttiker Tap(E) up 
to normalization A. We also verified that DOS taken from 
O atoms in other choices of (L,R) pairs yield agreement 
between Tap(E) and T^ P (E) comparable to that shown 
in Fig. EJa). Thus, the trustworthiness of Eq. (TTJ is well 
established, and we can apply it with justification to the 
FP-LSDA and QSGW Hamiltonians. 

The right panel of Fig. shows T% p {E), assuming the 
same normalization A (0.20 Ry 2 ) obtained by matching 
T% P (E) to Tap(E) in the ASA-LSDA approximation. 
The ASA calculation from the left panel is redrawn (blue 
solid line), and compared to a full-potential LSD A result 
(green dashed line). More precisely, the ASA approxima- 
tion to DOS£ 1 £ n (0| ? ) is replaced by its analog calculated 
with a FP-LSDA method. Finally we obtain Tf P (FP- 
QSGW) from DOS^ n (0^) calculated by QSGW in the 
repeated-slab geometry with 5 Fe and 5 MgO layers. 

As Fig. [2] shows, the peak in T® P falls near V in 
both the ASA-LSDA and FP-LSDA cases, in a gree ment 
with earlier full-potential LSDA calculations [la, Il7 |. 
The peak of the QSGW-derived transmission is shifted 
to higher energy by approximately 0.12 eV relative to 
the Ta p (FP-LSDA) result, putting it in close correspon- 



dance with the peak in \d 2 I/dV 2 \ (measured in APC, 
shown as black dashed line on Fig. [2j \d 2 I/dV 2 \ is shown 
rather than dl/dV, since the features are more easily seen 



12|.) Note that T% P (Q,SGW)<T% P (FP-LSDA). This is 



because the LDA gap (4.7 eV) is much smaller than the 
experimental (7.8 eV) and QSGW (8.8 eV) Q gaps; con- 
sequently 7k_E is overestimated in the LDA. 

Since finite-size effects of the leads can be impor- 
tant, we checked their sensitivity by calculating T^ P (FP- 
LSDA) for slabs with 5, 7, and 9 Fe layers and 5, 7, 
and 9 MgO layers. We found that the peak position and 
general shape of T® p depends weakly on the number of 
Fe and MgO layers. This is because minority IRS are 
mostly localized near the interface and their DOS quickly 
converges with number of Fe and MgO layers. On the 
other hand, finite size effects alter the majority DOS in 
the repeated-slab geometry. To eliminate finite-size effect 
in majority DOS we used the TB-LMTO-ASA majority 
DOSj^ J (0|') obtained for the semi-infinite electrode ge- 
ometry and thick, 12-layer, MgO barrier for all three cal- 
culations. Since DOS^^O^) is almost independent of 
energy on the scale we consider here (Ep±0A eV), the 
shape and peak position of the T^ p (E) will not depend 
on whether ASA-LSDA, FP-LSDA, or FP-QSGWis used 
to evaluate DOS^ J (03) in the semi-infinite limit. 

Significantly, there is only one peak in T in both the 
LSDA and QSGW approximations, despite the fact that 
two distinct peaks in the minority DOS appear (Fig. [TJ. 
To explain why only a single peak is seen, we analyze 
the QSGW DOS resolved by k in the 2D Brillouin zone 
of the (001) plane. Fig. [3] shows the k resolved DOS, at 
Ep + 0.12 eV and E F + 0.26 eV (see two peaks in Fig. [TJ. 
As Fig. [3] shows, the k-resolved DOS at E F + 0.12 eV 
is located mainly around the T point (k=0), while the 
DOS at the higher-energy peak {Ep + 0.26 eV) is concen- 
trated near the zone boundary (this feature is common 
to QSGW and LDA). The conduction band of MgO is 
free-electron like, thus imaginary wave number depends 
on k approximately as 7(fc) « y/k 2 , + k 2 [13j. Nearly all 
the DOS weight for the high-energy peak occurs at large 
k where 7 is large, so its contribution to T is effectively 
extinguished. This is in contrast to the low energy peak, 
where the surface DOS is concentrated at small k. 

In conclusion, we developed a parametrization that re- 
liably captures the essential elements of the tunneling 
process in magnetic tunnel junctions for any one-body 
hamiltonian. The tunneling probability derived from the 
QSGW approximation are in excellent agreement with 
observed differential conductance, in contrast to LSDA 
results. This confirms that the measured differential 
conductance peak is an intrinsic property of the ideal 
Fe/MgO (001) interface. This work also shows that corre- 
lations treated in the QSGW approximation are sufficient 
to realistically describe such metal insulator interfaces. 
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FIG. 2: (color online). Transmission function in the Fe/MgO system. Left: T AP (E), calculated in the APC by the TB-LMTO- 
ASA method for Fe/MgO/Fe MTJ with four MgO layers (red solid line) compared to the analytic formula T°p{E), Eq. JT} (blue 
dashed line), for a particular choice of sites in the MgO where local DOS is calculated. The latter was scaled by A=0.20Ry 2 . The 
quality of agreement is a measure of the quality of approximations used to obtain Eq. Jl}, as described in the text. We verified 
that the close correspondence between the analytic formula and the Landauer-Buttiker transmission seen in the left panel is 
insensitive to the choice of site chosen for the local DOS. Right: T AP (E), Eq. {l}, calculated in the antiparallel configuration for 
several cases. T®p shown in the left panel for the TB-LMTO-ASA method is redrawn as a (blue) solid line. The full-potential 
LSDA result is shown as a (green) dashed line. ^(Blfrom QSGW is shown as a (red) dash-dotted line. Experimental data 
(black dotted line) shows \d 2 I/dV 2 \ measured in APC [12] (arbitrary units). 




FIG. 3: (color online). The k-resolved DOS of minority elec- 
trons projected to surface Fe layer calculated by QSGW for 
energies corresponding to the peaks in QSGW DOS, Fig. 1. 
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